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Abstract. Wc introduce a new machinery to study the large time behavior for 
general classes of Hamilton- Jacobi type equations, which include degenerate parabolic 
equations and weakly coupled systems. We establish the convergence results by using 
the nonlinear adjoint method and identifying new long time averaging effects. These 
methods are robust and can easily be adapted to study the large time behavior of 
related problems. 



1. Introduction 

In this paper we introduce a new method to study the large time behavior for gen- 
eral classes of Hamilton-Jacobi type equations, which include degenerate quasilinear 
parabolic equations and weakly coupled systems. We consider nonlinear evolution 
Hamilton-Jacobi equations of the following abstract form 



ut + n{z,[u\) =0, int/x(0,oo) 
u{z, 0) = uo{z), on U. 



In order to present a general framework, including also the case of systems, we use 
a compact notation as follows. The set U is defined as f/ = T" x {1, 2, . . . , A;} with 
n,k & N fixed. The initial condition uq : U — )■ M and the Hamiltonian H : U x 
M X M" X M"y^ — )■ M are given, while u : U x [0, oo) -> M is the unknown. Here 
^syrn dcuotes the set of symmetric n x n real- valued matrices, and we set [f]{z,t) = 
{f{z,t),D^f{z,t),Dl^f{z,t)) for every smooth function / : f/ x [0, cx)) — )■ M. When 
needed, for every z ^ U we will also use the decomposition z = {x,i), with x G 
and 2 G {1,2,..., k}. We are only dealing with viscosity solutions of Hamilton-Jacobi 
equations and thus the term "viscosity" will be omitted henceforth. 

Our goal is to establish a unified method to prove that, under appropriate conditions, 
the following convergence result holds true 

max \\u{-,i,t) — {v{-,i) — ct)\\ioo(jn^ ^ as t — )■ oo, (1.1) 

i=l,...,k 
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where the pair [v, c) G C(T")^ x M is a solution of the ergodic problem 

{£) H{z, [v]) = c, in U. 

Under fairly general conditions, which will be specified later, there exists a unique 
constant c € M such that (S) holds true for some v G C{T^)^. It is however typically 
the case that such v is not unique, even up to additive constants, which makes the 
convergence (11 .ip delicate to be achieved. We give two important applications in 
this paper, one of which includes an open problem in the study of the large-time 
asymptotics for Hamilton-Jacobi type equations. More precisely, we study degenerate 
viscous Hamilton-Jacobi equations in Section 2 and the weakly coupled systems of 
Hamilton-Jacobi equations in Section 3. 

Before giving more details about the main ideas of the paper, let us briefly recall the 
state of the art of the subject. In the last decade, the large time behavior of solutions 
of single Hamilton-Jacobi equations, 

ut + H{x,Du) = 0, inrx(0,oo), (1.2) 

where H is coercive, has received much attention and general convergence results for 
solutions have been established. The first general result in this direction was proven 
by Namah and Roquejoffre in [15], under the additional assumptions: p t— )■ H{x,p) is 
convex, and H{x,p) > H{x,0) for all {x,p) G T" x M" and meiXxfzjn H{x,0) = 0. The 
same type of convergence results was obtained by Fathi [9] using a dynamical system 
approach from weak KAM theory. In this case, smoothness and uniform convexity on 
H{x,-) are assumed, i.e., DppH{x,p) > al for all (x,p) G x M" and a > 0, but 
no specific structure like the one above of Hamiltonians is needed. Afterwards, Davini 
and Siconolfi in [6] and Ishii in [10] refined and generalized the approach of Fathi, and 
studied the asymptotic problem for Hamilton-Jacobi equations on T" and the whole n- 
dimensional Euclidean space, respectively. Besides, Barles and Souganidis [2] obtained 
this type of results, for possibly non-convex Hamiltonians, by using a PDE method in 
the context of viscosity solutions. Barles, Ishii and Mitake [1] simplified the ideas in 
[2] and presented the most general assumptions (up to now) to achieve (11.11) . 

Very recently, Mitake and Tran [13] have been able to establish large time behavior 
results for weakly coupled systems of Hamilton-Jacobi equations in the general setting, 
by using a new representation formula for the solutions and the dynamical approach, 
which was inspired by the papers [6], [lO] . Also we refer to [I2] for related studies. After 
the work [13], Nguyen [13] generalized PDE techniques in [2] to obtain similar results. 

For the viscous Hamilton-Jacobi equation, 

ut + H{x,Du) = Au, inT"x(0,oo), (1.3) 

Barles and Souganidis [3] got the same type of asymptotics (II. ip . In this case, the 
associated ergodic problem has a structure which is much simpler than that one of 
(11.21) . Indeed, the strong comparison principle holds, and therefore the solution v of 
{£) is unique up to constants. The proof of (11. ip in [3] actually strongly depends on 
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this fact. We emphasize that though the large time behaviors are of the same type, the 
methods of obtaining convergence results for ( 11. 2p and ( 11. 3p are completely different. 

This paper presents a new method of obtaining large time behavior (11. ip of (C) based 
on the unified principle, which is to prove in the viscosity sense that 

\\ut{-,ht)\\L°°{T^) as t — > oo, for alH G {1, . . . , /c}. (1.4) 

Indeed, the stability of viscosity solutions, together with (11. 4p . suffices to achieve the 
convergence (11.11) . We notice that this principle already appeared in the papers of 
Fathi [9], Barles and Souganidis [2], in an implicit way. 

Let us briefly describe the key ideas of this paper. Without loss of generality, we 
assume the ergodic constant c = henceforth. In order to understand the limit as 
t — >■ oo, we introduce a rescaled problem. For e > 0, set u^{z,t) = u{z,t/e). Then 
{u'^)t{z,t) = e~^{u)t{z,t/e), Du'^{z,t) = Du{z,t/e), and solves 



euf + n{z,[u^]) = 0, in?7x(0,oo) 
u^{z,0) = Uo{z), on U. 



By this rescaling, u'^{z,l) = u{z,l/e) and we can easily see that to prove ( 11. ip is 
equivalent to prove that 

\\u'^{-,i,l) ^ v{-,i)\\L^(jn) as e — )■ 0, for alH G {1, . . . , A;}. (1.5) 

To show (11. 5p . we first introduce the following approximation of (C)^ 

ewf + n{z,[w^]) = e^Aw^, inf/x(0,oo), 
w'^{z,0) = Uq{z), on U. 

Then, we observe that w"^ is smooth, and 

||ii;'^(-, z, 1) — u^(-, z, l)||ioo(T") — )■ as e — )■ 0, for alH G {1, . . . , /c}. 

It is thus enough to derive the convergence, as 5 — 0, of w^(-, i, 1) by proving that 

e\\wl{-,i,l)\\L°^(jn'j as e — 0, for alH G {1, . . . , /c}, (1.6) 

which is actually a rigorous version of (11. 4p . 

In order to prove (11. 6p . we combine the nonlinear adjoint method introduced by 
Evans [7] and a new observation on the averaging action as t — )■ oo (or equivalently as 
e — )■ by rescaling), as clarified below. Let be the formal linearized operator of 
{A)£ around w^, i.e. 



C^^f := ^ {e{w' + Vf)t + K + Vf])- ^'A(u;^ + r//) 



»7=o' 



for any f E C^{U x (0, oo)). We consider the following adjoint equation 



Cl^a' = 0, in f/ X (0, 1) 
aHz, I) = 5,,„, on 



4 



F. CAGNETTI, D. GOMES, H. MITAKE, H. V. TRAN 



where C*^e is the formal adjoint operator of C^^, and 5^^ is the Dirac delta measure at 
some point zq G U. We then see that cr^(-,t) is a probability measure for all t G (0, 1) 
and conservation of energy holds, namely 

dt Jjn 

The conservation of energy in particular gives us a different way to interpret ewf{-, 1) 
as 

pi p ^ 

ewl{zo,l)= / ^ K]) -e^Au;")a"(x,z,t)rfxdt. (1.7) 

Jo Jt" 

We then use techniques of the nonlinear adjoint method to establish that the right 
hand side of (11. 7p vanishes as £ — ?■ 0, which is basically equivalent to obtain the main 
principle (11.61) . Notice that the averaging action appears implicitly in (11.71) and plays 
the key role here. More precisely, if we rescale the above integral back to its actual 
scale, it turns out to be 

[ V (?/(x,z, [M;^](a;,£t)) -e^Atf;^(x,z,£t))a^(x,z,£t)cixc?t, (1.8) 
T Jo Jt" j^j^ 

where T = 1/e — )■ oo. 

The nonlinear adjoint method for Hamilton- Jacobi equations was introduced by 
Evans [7] to study the vanishing viscosity process, and gradient shock structures of vis- 
cosity solutions of non convex Hamilton-Jacobi equations. Afterwards, Tran [17] used 
the nonlinear adjoint method to establish a rate of convergence for static Hamilton- 
Jacobi equations and was able to relax the convexity assumption of the Hamiltonians 
in some cases. Cagnetti, Gomes and Tran [1] then used it to study the Aubry-Mather 
theory in the non convex settings and established the existence of Mather measures. 
See also [H |8] for further developments of this new method in the context of Hamilton- 
Jacobi equations. 

What is new in this paper is an observation based on the combination of the nonlinear 
adjoint method with the averaging action as t -> oo presented roughly above, which 
gives a new method to achieve the convergence (11.11) . This gives not only a unified 
approach on the study of large-time asymptotics for (II. 2p . (II. 3p and weakly coupled 
systems but also an answer for that of degenerate viscous Hamilton-Jacobi equations 
with convex Hamiltonians. Since our purpose is largely expository, we have chosen 
only two important applications which we have already mentioned, but the method is 
robust and can easily be adapted to study the large time behavior of related problems. 
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2. Degenerate Visgous Hamilton- Jagobi equations 

In this section we study degenerate viscous Hamilton-Jacobi equations. To keep the 
formulation as simple as possible, we consider an equation of the form 

(ut + H{x,Du) = a{x)Au, in T'^ x (0, oo), 
\u{x,0) = uo{x), on T*^, 

where Mq G C(T"). Throughout this section we always assume the following. 
(HI) He C2(T" X R"), and there exists e>0 such that 

(H2) There exists C > such that 

\D,H{x,p)\<C{l + \p\'). 

(H3) a G C2(T"), a > 0. 
We notice that if a = or a = 1, then the equation in (C) becomes the first-order 
Hamilton-Jacobi equation (11.21) . or the viscous Hamilton-Jacobi equation (II. 3p . re- 
spectively. We emphasize that in general, we can deal with possibly degenerate viscous 
Hamilton-Jacobi equations at some part, i.e., {a = 0} is not empty and can be arbi- 
trary. The next three propositions concern basic existence results, both for (C) and for 
the associated stationary problem. The proofs are standard, hence omitted. 

Proposition 2.1. Let uq G C(T"). There exists a unique solution u of (C) which 
is uniformly continuous on T" x [0,oo). Furthermore, if Uq G Lip(T"), then u G 
Lip(T" X [0,cx))). 

We next state the basic existence result on the stationary problem. 

Proposition 2.2. For every e G (0, 1) there exists a unique constant such that the 
following ergodic problem 

H{x, Dtf) = {e^ + a{x))Av' + 11^, in T" (E)^ 

has a unique solution if G Lip (T") up to constants. In addition, 

[He -c\< Ce\ \\Dv'\\Lo.(f^) < C, (2.1) 

for some c G M and a positive constant C independent of e. 

In view of the quadratic or superquadratic growth of the Hamiltonian if, we can 
get (12. ip by the Bernstein method. By passing to some subsequences if necessary, 

— v^{xq) for a fixed xq G T" converges uniformly to a Lipschitz function f : T" — M 
which is a solution of 

H{x, Dv{x)) = a{x)Av + c, in T". (E) 

Proposition 2.3. There exists a unique constant c G M such that the ergodic problem 
(E) admits a solution v G Lip (T*^). 
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Without loss of generality, we may assume c = henceforth. We can now state the 
long time convergence result for (C). 

Theorem 2.4 (Main Theorem 1). Letu be the solution of (C) withuo € C(T"). There 
exists a solution v of (E) such that 

t) — t'||Loo(X") — as t +00. 

The proof of Theorem 12.41 is the subject of the following subsections. 

2.1. Regularizing Process and Proof of Theorem 12. 4L Without loss of generality 
we can assume that uq G Lip (T"), since the general case uq € C(T") can be obtained 
by a standard approximation argument. In particular, thanks to Proposition 12.11 we 
see that u is Lipschitz continuous on T" x [0, oo). 

As stated in Introduction, we consider a rescaled problem. Setting u'^{x,t) = 
u{x, t/e) for e > 0, where u is the solution of (C), one can easily check that satisfies 



(C)e 



ul + H{x,Du^) = a{x)Au\ in x (0, oo) 
0) = no(x), on T". 



Notice however that in this way we do not have a priori uniform Lipschitz estimates 
on e, since the Lipschitz bounds on u give us that 

ll'"tlU°°{T"x[o,i]) < C*/^; ||-Dn^||ioo(Trnx[o,i]) < C. (2.2) 

In general, the function is only Lipschitz continuous. For this reason, we add a 
viscosity term to (C)^, and we consider the regularized equation 



ewt'"^ + H{x, Dw^^'i) = (a(x) + r])Aw^''^, in T" x (0, oo), 
w^'''(x, 0) = uo{x), on T" 



for r] > 0. The advantage of considering (A)^ lies in the fact that the solution is smooth, 
and this will allow us to use the nonlinear adjoint method. We specially denote (A)^ 
by (A)e, and w^'^ by w'^. 

The proof of Theorem 12.41 is based on the two following results. 

Lemma 2.5. Letw'^ he the solution of (A)g. There exists C > independent of e such 
that 

lk'(-, l)||cHT") < C, \\u%-, 1) - l)|Uc.(Trn) < Ce. 

The above Lemma is a basic result on the rate of convergence of the regulariz- 
ing process. In principle, the rate of convergence for the vanishing viscosity method 
is \/ viscosity coefficient /(the coefficient of wf) and therefore, we naturally expect the 
above result. The proof of Lemma [2.51 can be derived using standard arguments. Nev- 
ertheless, we sketch the proof below, since some of the estimates involved will be used 
later. 

Before starting the proof, we state a crucial property of the function a G C'^iT^). 
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Lemma 2.6. There exists a constant C > such that 

\Da{x)\'^ < Ca{x) for all x E T". (2.3) 



Proof. In view of [161 Theorem 5.2.3], a^^^ G Lip (T"). It is then immediate to get (12. 3 p 
by noticing that, for a{x) > 0, 

Proof of Lemma 12.51 Since we are assuming that the ergodic constant for (E) is now 
and H is quadratic or superquadratic on p, we can easily get the first estimate. We 
only prove the second one. 

Let w^''' be the solution of (A)^ and set 4>{x,t) = \Dw^''^\'^ /2. Then, satisfies 

e(t)t + DpH ■ D(f) + D^H ■ Dw''"^ = ij] + - (r/ + a)\D'^w'^'^\^ + a^.^^^^Aw"'''. 



We next notice that for 5 > small enough, in light of Lemma 12.61 

a^^wi'^^Aw''"^ < C\Da\ ■ |Au;"'''| < ^ + 5\Da?\D^w''''\^ <C + -a\D^w'''^\'^ . 
Hence 

e(t)t + DpH ■ D4> + ^{r] + a)\D^w'''^\^ < (r/ + a)A0 + C. (2.4) 
Multiply (12.41) by and integrate over [0, 1] x T" to yield 

[ [ {a{x)+r])\D^w''''\^a'dxdt<C (2.5) 

for some C > 0. 

Note that w'^'^ is differentiable with respect to r] by standard regularity results for 
elliptic equations. Differentiating the equation in (A)^ with respect to rj, we get 

e (u;^'")^ + DpH{x, Dw''"^) ■ Dw'^ "^ = Aw''" + {a{x) + r/) Au;^'", in T", 

where denotes the derivative of a function / with respect to the paramenter 77. 
Multiplying the above by a'^'" and integrating by parts yield 

ew'^'\xo,l)= [ [ Aw^'V^'^rfxrft, 

where we used the fact that w^'"{x, 0) = 0. Thanks to (12. 5p . by the Holder inequality 

e\w'^''^{xo,l)\<c([ [ I^V'^lV^'^rfxcit) < 



By choosing properly the point xq we have thus 

IK(-,i)IU^(T")<^, 

which gives 

||w;^'^(-,l)-n^(-,l)|Uo.(Tn)<^l^. 
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Finally, observing that = w^'^*, choosing t] = we get the result. 

□ 

Next theorem gives (11.61) in the special case of problem (A)^. 
Theorem 2.7. We have 

\ime\\wl{-, l)llLoo(Trn) = 0. 

More precisely, there exists a positive constant C, independent of e, such that 

e\\wt{-, l)IUoo(Trn) = \\H{; Dw%-, 1)) - {e" + a{-))Aw'{; < Ce'/\ 

The proof of Theorem 12.71 is postponed to the end of this section. We can now give 
the proof of Theorem 12.41 

Proof of Theorem \2A[ By Lemma [2.5[ we can choose a sequence {Sm} — ^ such that 
{w'^"^{-,l)} converges uniformly to a continuous function v. By Theorem 12. 7[ w is a 
solution of (E), and thus a (time independent) solution of the equation in (C)^. We let 
tm = l/£m and use Lemma 1231 to deduce that 

tm) — 'y||L°°(T") — > as m — oo. 

Let us show that the limit does not depend on the sequence {tm}mgN- Now, for any 
a; G T", t > such that t^ <t < tm+i, we use the comparison principle to yield that 

\u{x,t) -V{x)\ < \\u{-,tm + (t-tm)) - ^K") IU°°(T") < ||M(-,^m) " ^^l") IU°°(T") • 

Thus, 

lim |-u(x, t) — f (x)| < lim ||m(-, tm) — 'i'(')IU°°(ir") = 0, 

t—^oo m— >oo 

which gives the conclusion. □ 

2.2. Convergence Mechanisms: Degenerate Equations. The adjoint equation 
corresponding to (A)^ is 



(AJ), 



-eaf - diY{DpH{x, Dw')a') = e^Aa' + A{a{x)a'), in T" x (0, oo) 
_a^(x, 1) = 5^0, on T", 

where is the Dirac delta measure at some point xq G T". 
Lemma 2.8 (Elementary Properties of a^). We have > and 



/ a''{x,t)dx = 1, for all t e [0,1]. 



This is a straightforward result of adjoint operator and easy to check. We show 
now the following key lemma, which provides integral bounds on first and second order 
derivatives of the difference — on the support of and a. 
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Lemma 2.9 (Key Estimates). There exists a positive constant C, independent of e, 
such that the following hold: 

\D(w' - + e'^lDHw' - v')\A a'dxdt < C, 



« /7 (ji 

Jo ^Jt" \t 

(ii) [ I a^{x)\D''{w' -v')\''dxdt<C^/i. 
Jo Jt" 



Proof. Subtracting equation (A)^ from (E)^, thanks to the uniform convexity of H, we 
get 

= e{v' - w% + H{x, Dv') - H{x, Dw') - {e^ + a{x))A{v' - w') - He 
> e{v' - w')t + DpH{x,Dw') ■ D{y' - w') + e\D{if - w")p 
- {e^ + a{x))A{v' - w') -H,. 

Multiply the above inequality by and integrate by parts on [0, 1] x to deduce 
that 

t=i 



e 



Diw' -v')\^(j' dxdt<He + e 



-'T" 



[w^ — tf)a^ dx 



t=0 



He + e{w^{xQ,l) — v'^{xq)) — e I {uq{x) — v'^{x))a^{x,Q) dx 

He + ew^{xQ,l) — £ / {tf{xQ)—v^{x))a^{x,Q)dx — £ / uo{x)a^{x,Qi) dx 



T" 



<He + Ce + CeWDv'^WL^ijn) — e I uo{x)a'^{x,0) dx < Ce, 

where w is a solution of (E) and we used Lemma 12.51 and Proposition 12.21 (recall that 
we set c = 0) in the last two inequalities. We hence get 



[ [ \D{w' -v')\^a'dxdt <Ce, 
Jo J J" 



(2.6) 



which is the first part of (i). 

Next, subtract (A)^ from (E)^ and differentiate with respect to Xi to get 



e{v' - w')^^t + DpH{x, Dif) ■ Dvl^ - DpH{x, Dw') ■ Dwl^ 

+ H,Xx,Dv') - H^X^, Dw') - {e^ + a)/\{v' - w'),^ - a,, A(t;^ - w') = 0. 

Let (p{x,t) = \D{v'^ — w'^)\'^ /2. Multiplying the last identity by {v'^ — w^)^^ and summing 
up with respect to i, we achieve that 



Eipt + DpH{x, Dw^) ■ Dip 



D„H(x. Dv') - D^Hix. Dw') ■ Dv 



Xi 



+ [d^H{x, Dv') - D^H{x, Dw')^ ■ D{y' - w') - {e^ + a(x))A^ 

+ {e^ + a{x))\D^{v' - w')\'^ - [Da ■ D{v' - w')] A{v' - w') = 0. 
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By using various bounds on the above, we derive that 

e^t + DpH{x, Dw') ■ Dip - {s^ + a{x))Aip + {e^ + a{x)/2)\D\v' - w") ^ 

<C{\Dh!'\ + l)\D{v' -w')\^ (2.7) 

The term in the right hand side of (12 .7^ is a dangerous term. We now take advantage of 
( 12. 5p and ( 12. 6p to handle it. Using the fact that ||Dt;^||ioo and ||Z^tf^||Loo are bounded, 
we have 

C|DV| \D{v'-w')\^< C\D\v'-w')\ \D{v' -w')\^ + C\D^w'\ \D{v' - w')\^ 

Combine (ET]) and <^Bj to deduce 
e^t + DpH{x, Dw') ■ D<p - {e^ + a{x))A<p + —\D\v' - w')\^ 



< -\D''{tf - w')\^ + ^\D{v' - w')\^ + C\D^w'\. (2.8) 



< C\D{v' - w')\^ + —^\D{tf - w;")|2 + C|DV|. (2.9) 



We multiply ( 12. 9p by a^, integrate over [0, 1] x T", to yield that, in light of ( 12. 5p and 



JT" ^ 



/" I \D^w'\a'dxdt 
Jo Jt" 



<£^ + c( [' [ \D^w'\^a'dxdt)"\ f f a'dxdtV^' <^ + ^<^ 
Jo JT" ^ ^ Jo Jt" 



£^ ^ JO J J" ' Wo T" ^ ^'^ 

Finally, we prove (ii). Setting ip{x,t) = a{x)\D{v'^ — w^)(x, t)p/2 = a{x)ip{x,t) and 
multiplying (12. 7p by a(x) we get 

eipt + /^pi/(x, Dw') ■ Dip - {DpH{x, Dw') ■ Da)ip - {e^ + a{x))Aip 

+ {e^ + a(x))(Aa(/5 + 2Da ■ Dip) + a{x){e^ + a{x) /2)\D^{v' - w')\'^ 

< Ca{x){\D^v^ + l)\D{v' - w')\\ 
We use the facts that Da, Aa are bounded on T" to simplify the above as follows 

ei/jt + DpH{x, Dw') ■ Dtp - [e^ + a{x))Ai^ + a{x){e^ + a[x)l2)\D'^{v' - w^)"^ 

< C\D{v' - w')\^ - 2{e^ + a{x))Da ■ Dip + Ca{x)\D\'\ \D{v' - w')\^. (2.10) 

Next, we have to control the last two terms on the right hand side of (I2.10p . Observe 
first that for 6 > small enough 

2\{e^ + a{x))Da ■ Dip\ < C{e^ + a{x))\Da\ \D\v' - w')\ \D{v' - w')\ 

< 6{e^ + a{x))\Da\^\D\v' - w')\^ + ^\D{v' - w')\^ 



< lie^ + a{x))a{x)\D^{v' - w')\^ + C\D{v' - w')\\ (2.11) 
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where we used Lemma 12.61 in the last inequahty. On the other hand, 

a{x)\D^v^ \D{if-w')\^ 
<a{x)\D^w'\ \D{v' - w')\^ + a{x)\D\v' - w')\ \D{v'-w')\ 

< ^/Ea{x)\D^w'\^ + -^\D{v' - w')\^ + ^^\D\v' - w')\^ + C\D{v' - w')\\ 

\ E 8 



(2.12) 

We combine (ElO]), f l2TT]l . and fl2:T2D to deduce that 

e^\)t + DpH{x, Dw') ■ - {e^ + a(x)) At^ + ^!^\D'^{v^ - w") P 

< (C + Ce-^'^)\D{v' - w')\'' + e^'^a{x)\D^w'\^. 

We muhiply the above inequahty by a^, integrate over T" x [0, 1] and use (12.51) to yield 
the result. □ 

2.3. Averaging Action and Proof of Theorem 12.71 

Lemma 2.10 (Conservation of Energy). The following hold: 
(i) / {H{x,Dw')-{e^ + a{x))Aw')a'dx = Q, 

dt Jjn 



ill) 5<(Xo,l) 



f [ {H{x,Dw')-{e^ + a{x))Aw')a'dxdt. 
Jo Jt" 



We stress the fact that identity Lemma [2. 101 (ii) is extremely important. As stated in 
the Introduction, if we scale back the time, the integral in the right hand side becomes 
(11.81) . that is the averaging action as t — >■ oo. Relation (ii) together with Lemma [2.91 
allow us to conclude the proof of Theorem 12.71 

Proof. We only need to prove (i) as (ii) follows directly from (i). This is a straightfor- 
ward result of adjoint operators and comes from a direct calculation: 

^ [ iH{x,Dw')-ie^ + a{x))Aw')a'dx 

dt Jjn 

= [ {DpH{x, Dw') ■ Dwl - {e^ + a{x))Awl)a' dx 

+ I {H{x,Dw') - {e^ + a{x))Aw')atdx 
= - (div {DpH{x, Dw')a') + (e^ + a{x))Aa')wl dx - ewlal dx = 0. □ 

We now can give the proof of Theorem 12.71 which is the main principle to achieve 
large time asymptotics, by using the averaging action above and the key estimates in 
Lemma 12. 9[ 
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Proof of Theorem 12.71 Let us first choose Xq such that 

\ew^{xo,l)\ = \H{xo,Dw%Xo,l)) - {e^ + a{xo))Aw%Xo,l)\ 

= \\H{; DW%; 1)) - is' + a{-))Aw'{; l)IUoo(Tn). 

Thanks to Lemma 12.101 and Proposition 12.21 

l)|U.o(Tn) = \\H{; Dw'i; 1)) - (s^ + a(-))Au;^(-, 1)|U-(T") 

= [ [ {H{x,Dw') - {e^ + a)Aw')a' dxdt 

Jo Jj'^ 

< [ [ \{H{x,Dw') - {e^ + a)Aw') - {H{x,Dv') - {e^ + a)Av')\a' dxdt+\H. 

Jo Jj" 

< [ [ [C\D{w' -v')\ + {e^ + a)\A{w' -v%a'dxdt + Ce^. 

Jo J J" 

We finally use the Holder inequahty and Lemma 12.91 to get 

<c([ [ \D{w' -v')\'^a'dxdt] +Ce^( [ [ \D\w' - v')\^a' dx dt) 

\Jo Jj" J \Jo JT" / 



rl r \V2 

+ C{ j / a'^{x)\D^{w' -v')\^(j'dxdt\ +Ce^<Ce^'\ □ 



3. Weakly Coupled Systems of Hamilton-Jacobi Equations 

In this section, we will basically repeat the steps which we described for the single 
equation in Section [2], adjusting to the weakly coupled system of Hamilton-Jacobi 
equations: 

' {ui)t + Hi{x, Dui) + ui-U2 = 0, in T'^ x (0, oo), 
(SC) {{u2)t + H2{x,Du2)+U2-ui = 0, inT"x(0,oo), 

Ui{x, 0) = Uoi{x), on T*^, for i = 1, 2. 

Throughout this section we always assume that u^i G C(T") and that the functions Hi 
satisfy (HI), (H2) for z = 1,2. 

We first state the basic existence results for (SC) and for the associated stationary 
problem. The proofs of the next three propositions are standard, hence omitted. 

Proposition 3.1. Let (^01,^02) G C(T")^. There exists a unique solution (^1,^2) of 
(SC) which is uniformly continuous on T" x [0,oo). Furthermore, if uqi G Lip(T'^), 
then Ui G Lip (T" x [0, 00)) for z = 1, 2. 
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Proposition 3.2. For every e > sufficiently small there exists a unique G M such 
that the following ergodic problem 



(SE), 



Hi{x, Dvf) + vf -vl = e^Avl + H^, in T", 
H2{x, Dvl) + vl-vl = e^^vl + H,, in T", 

has a unique solution G Lip (T")^ up to additional constants. In addition, 

\He-c\<Ce^, P<IU-(T") < C, « = 1,2, 

for some c G M and a positive constant C independent of e. 

The following result is straightforward. 

Proposition 3.3. There exists a unique constant c G M such that the ergodic problem 

Hi{x, Dvi) + f 1 — f 2 = c, in T", 
H2{x, DV2) + V2 - vi = c, in T", 



(SE) 



has a solution [vi, V2) G Lip (T 

Without loss of generality, we may assume c = henceforth. We present our second 
main result. 

Theorem 3.4 (Main Theorem 2). Let {ui,U2) be the solution of (SC) with (moi, Mo2) ^ 
C(T")^. There exist a solution {vi,V2) G (C(T"))^ such that 

\\ui{-,t) — ViWioo^jn) ^ as t — 7- +00, fori = 1,2. 

3.1. Regularizing Process and Proof of Theorem 13.41 Without loss of generality, 
in the following we will assume that {ui, U2) is Lipschitz on T" x [0, 00), as it was done 
in Section O Once again, we will follow the machinery stated in the Introduction. 

We perform a change of time scale. For e > 0, let us set ul{x,t) = Ui{x,t/e), which 
is the solution of 

e{ul)t + Hi{x, Du\) + u\-ul = Q, in T" x (0, 00) 
(SC), {e{ul)t + H2{x,Dul)+ul-u\ = Q, in T" x (0, cx)) 

uf (x, 0) = UQi{x), on T", for z = 1, 2, 

and we approximate (SC) by adding viscosity terms to the equations 

e{w\)t + Hi{x, Dw\) + wl-wl = e^^w\, in T" x (0, oo), 
(SA), {e{wl)t + H2{x,Dwl) +wl-wl = e^/^wl, in T" x (0, cx)), 

wl{x, 0) = UQi{x), on T*^, for i = 1,2. 

We can conclude the proof of Theorem 13.41 with the following two results. 

Lemma 3.5. Let {w{,w'^ be the solution of (SA),. There exists C > independent 
ofe such that ||u;f (-, 1)||ci(T") < C, ||uf(-,l) - w;f (-, 1) ||l-(t") < Ce fori = 1,2. 
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Theorem 3.6. We have 

limmax£||«)t(-, l)||Loo(Trn) = 0. 

£— 1=1,2 

3.2. Convergence Mechanisms: Weakly Coupled Systems. The adjoint system 
corresponding to (SA)^ is 

{-e{al)t + div{DpHi{x, Dwl)al) + - a'^ = £^Aaf , in T" x (0, oo), 
+ diY{DpH2{x, Dwl)al) + a'2 - af = e^Aaf, in T" x (0, 00), 
erf (x, 0) = Sik6:^o, on T", for z = 1, 2. 

where 6ik = 1 Hi = k and = if z 7^ k, and Xq G T*^ and A; G {1, 2} are to be chosen 
later. Notice that for any choice of k, either (■, 1) = or o'2{-, 1) = 0. Let us record 
some elementary properties of (af , a|) first. 

Lemma 3.7 (Elementary properties of {af, (t|)). We have cr| > for 2 = 1,2 and 

y2 / erf (x, t)dx=l for all t G [0, 1]. 

We next derive key integral bounds for (tyf , w^, (ff , V2) and their derivatives on the 
supports of (crfjCTg). 

Lemma 3.8 (Key estimates for weakly coupled systems). The following hold true: 

(ii) f [ [{wl-vl)-{w'2-v'2)]\al + a'2)dxdt<Ce. 
Jo Jt" 

Lemma 13.81 (ii) is a new observation on the study of weakly coupled systems, which 
gives us the large time average control on the coupling terms. This is actually the key 
point in the derivation of the main result for systems (Theorem 13.41 and Theorem 13. 6p 
as one can see in the proof of Lemma 13.91 

Proof. We will only prove (ii), since part (i) can be derived by repeating the proof of 
Lemma 12. 9[ 

Thanks to Lemma 13.51 we can always add to the pair (ff , v\) an arbitrarily large 
constant C (independent of e) such that 

2C > vl > wl in T", for 2 = 1, 2. (3.1) 

Let 0j = {vf — wfY /2 for i = 1,2. Subtract the first equation of (SA)^ from the first 
equation of (SE)^, and multiply the result by ff — wl to get 

e{vl - wl){vl - wl)t + {I'l - wl){H,{x,Dvl) - H,{x,Dwl)) 
+ {vl - wlf - {vl - wl){v'2 - wl) = e\vl - wl)A{vl - wl) + Hs{vl - wl). 
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We employ the convexity of Hi and (13. ip to deduce that 



e{(j)i)t + DpHi{x,Dwl) -04)1 



n — V2 



(3.2) 



Similarly, 



e{<p2)t + DpH2{x, Dwl) ■ D4>2 + ^2 - 4>i 



(3.3) 



Multiplying 03. 2p . f l3.3p by af, o"| respectively, and integrating by parts 



1 
2 



[{vl - wl) - {vl - w'^)fial + (t|) dx dt 



'0 JT' 
2 



bial dx 



+ C\He\-e^Y^ / \D{vl-wl)\^aldx<Ce, 



which implies (ii) 



□ 



3.3. Averaging Action and Proof of Theorem 13.61 For each i G {1,2}, setting 
j = 3 — z we have {i,j} = {1,2}. The following result concerning conservation of 
energy and averaging action is analogous to Lemma 12.101 and therefore we omit the 
proof. 



Lemma 3.9 (Conservation of Energy for weakly coupled systems). The following hold: 



d 
dt 



[ y^{Hi{x, Dwl) + wl - w'j - e^AwDal dx = 0. 



(ii) -{k- l)e{wl)t{xo, 1) - (2 - k)e{wl)t{xo, 1) 

2 



f f ^{Hi{x, Dwl) + wl - w'j - e^AwDal dx dt. 
Jo Jt" 



Proof of Theorem 13.61 Without loss of generality, we assume that there exists xq G 
such that 



£\iwl)tixo,l)\ = emax||«)t(-,l)||Loo(Tn). 

1=1,2 
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We then choose /c = 1 in (SAJ)^ and use Lemma [3^ to get 



e max 



/ / ^{H,{x,Dwl) + wl-w'^-e^Awl)aldxdt 

Jo Jjn 



< 



'- — Wj — e^Aw^) 



< 



ax||(u'-)t(-,l)IU-(T") 
r / y2\H,ix,Dw 

- {Hi{x, Dvl) + vl - if. - e'AvX)^Gl dx dt 
f [ J2[C\Di^^-^t)\+e'\A{wt-vt)\]aldxdt 

Jo Jt" -^-^ 



+ \He 



+ I [ \ [{wl - vl) - {w'^ - v'^)] K - \dxdt+\H, 
'0 Jt" 



Thus, 



/ / y2'^H,{x,Dwl)+wt-w'j-e^Aw',)a',dxdt 
Jo Jt" ^^-^ 



i=l 
2 



i=l 



1 

1 



\D(w^ -v'i)\'^a'idxdt 



JT 



1 




2 






\fh 




Jo Jt" 




1/2 



D'^iw', -vl)\^aldxdt 



\{wl-vl)-{wl-vl)\'aldxdt 



where the last inequahty follows by using Lemma 13.81 



□ 



3.4. Generalizations. It is also fairly straightforward to generalize the above system 
of two equations to systems of m-equations as follows. We consider a system of the 
form 

m 

{ui)t + Hi{x, Dui) + ^ djUj = 0, in T" x (0, 00). 
j=i 

and assume that Hi satisfy (HI), (H2) and ca > 0, Cjj < for i ^ j, Yl]jLi (^ij = 0- 

In order to achieve large time behavior results, we only need to establish the key 
estimate concerning the coupling terms as in part (ii) of Lemma 13.81 More precisely, 
we show that 

hm/ / J2\c,,\[{w^^-zf^)-{wt-vt)]'aldxdt = 0. (3.4) 
^^"Jo J J" 
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Set (pi = {vf — /2 for z = 1, . . . , m. Then we can compute that 
e{vl - wl){vl - wt)t + « - wt)m^. Dvl) - DwD) 

m 

The last term in the right hand side of the above identity can be written as 

m 

^Ci,{vl - wl){v'^ - w]) = \c,j\ {{vl - wlf - {vl - wl){v'^ - w'^^ 

= E ^^^■'^i + 2 5^ ''''^■l [^""j' ~ "^^^ ~ ~ 
i=i j=i 

Hence 

m 

+ DpHi{x, Dwl) ■ Dcpi + E Cij'/'j 

^ m 

Then ( 13 ■4p follows immediately. 

Finally, it is worth to mention that if we combine the arguments in Sections 2 and 3, 
then we can immediately get the result on the large-time behavior for weakly coupled 
systems of degenerate viscous Hamilton-Jacobi equations. We omit the details here. 
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